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ABSTRACT 

We expand the gauge field in terms of a suitably constructed complete set of 
Bloch wave functions, each labeled by a band designation n and a wave number 
K restricted to the Brillouin zone. A noncompact formulation of lattice QCD 
(or QED) can be derived by restricting the expansion only to the O^'^-band 
{n = 0) functions, which are simple continuum interpolations of discrete values 
associated with sites or links on a lattice. The exact continuum theory can be 
reached through the inclusion of all n = and n ^ bands, without requiring 
the lattice size i . This makes it possible, at a nonzero i , for the lattice 
coupling ge to act as the renormalized continuum coupling. All physical results 
in the continuum are, of course, independent of i . 
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1. INTRODUCTION 



In a recent paper-*^ (hereafter referred to as I) we proposed a new lattice 
formulation of the continuum field theory. Instead of the usual Fourier series, the 
continuum field operators are expanded in terms of a suitably chosen complete 
set of orthonormal Bloch functions 

{fn{K\^} (1.1) 

— * 

where K denotes the Bloch wave number restricted to the Brillouin zone, and 
n labels the different bands, similar to the one-particle wave functions in a 
crystal (i.e., e~'^^''^ fn{K \f) has the periodicity of the lattice). The lattice 
approximation is then derived by either restricting it to only one band (say, n = 
0), or to a few appropriately defined low-lying bands. Since the inclusion of all 
bands is the original continuum problem, there is a natural connection between 
the lattice and the continuum in this method. By including the contributions due 
to more and more bands, one can systematically arrive at the exact continuum 
solution from the lattice approximation. There is a large degree of freedom 
in choosing the Bloch functions (1.1), as the original continuum theory has no 
crystal structure. These extra degrees of freedom are analogous to gauge fixing; 
the final answer to the continuum problem is independent of the particular choice 
of Bloch functions (or the lattice structure). 

In paper I we gave a few simple examples of such functions and applied 
them to spin-0 and spin-i fields. We noted that this approach provides an 
effective means of removing the spurious lattice fermion solutions for the Dirac 
equation. The underlying reason is quite simple: After all, in any crystal in 
nature all the electrons do move in a lattice and satisfy the Dirac equation; yet 
there is not a single physical result that has ever been entangled with a spurious 
fermion solution. Consequently, it should not be that difficult to get rid of these 
unphysical elements. 

On a discrete lattice, particles hop from point to point, whereas in a real 
crystal the lattice structure is embedded in a continuum and electrons move 
continuously from lattice cell to lattice cell. In a discrete system, the lattice 
functions are defined only on individual points (or links, as in the case of the 
gauge field). However, in a crystal, each Bloch wave fn{K\f) is a continu- 
ous function in r , and herein lies one of the essential differences. In the case 
of gauge theories, there is an additional complication, since the usual lattice 
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gauge action^ is compact, which makes it intrinsically different from the non- 
compact continuum action (except in the limit of weak coupling or infinitesimal 
lattice size). As we shall see, these differences can be resolved by extending 
the lattice formulation developed in I to both quantum electrodynamics (QED) 
and quantum chromodynamics (QCD). The result is a new noncompact lattice 
formulation of gauge theories that can serve as the first approximation of the 
continuum theory in a systematic way, without requiring the lattice size £ — > . 

The freedom in choosing any complete set of Bloch functions allows us 
to explore a family of new functions that are convenient to use. For higher 
bands, we would like the Bloch functions to stay very close to the Fourier series 
of high wave numbers, so that renormalization calculations can be carried out 
analytically. On the other hand, the zeroth band functions fo{K \ f) should be 
made of simple continuum interpolations of the discrete values given at lattice 
(or link) sites; this way, the restriction to zeroth band would naturally resemble 
the usual lattice theory. Furthermore, in order to retain the "locality" character 
of the original continuum theory in the one-band approximation, we require 
fo{K If) to satisfy, for any lattice site j located at fj , 

1 ^e-^^-^^/o(i?|f) = ff \f-rj\ > vi (1.2) 

K 

with V a finite number 0(1) , i the unit lattice size, O the total volume » 
and K being summed over all wave numbers within the Brillouin zone. This 
lattice locality condition rules out the simple Fourier series for /o(i? | r ) , since 
for /o(-^ I r) = e*^'^ the corresponding sum (1.2) for a cubic lattice with 
ri = would be 

sin j{x- Xj) sin j {y - Vj) sin j {z - Zj) 
L3 sin l{x- Xj) sin ^{y- Vj) sin f (z - Zj) 

which is non-local. Instead, condition (1.2) leads to the concept of "lump func- 
tions" and "link functions", which will be introduced in Section 2. 

For applications to gauge theories, there is also the problem of compatibility 
between the continuum gauge-fixing condition and the band decomposition. We 
discuss QED first in the time-axial gauge in Section 3.1 and then in the Coulomb 
gauge in Section 3.2. By using suitable combinations of these lump and link 
functions we show that in both cases the continuum gauge requirement. 
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A4f,t) = (1.4) 

for the time-axial gauge or 

V-A{f,t) = (1.5) 

for the Coulomb gauge, can hold within each n^^ band at all r and t , where 
Afj_{r',t) is the electromagnetic field operator. Thus, as we shall see, in the 
lattice approximation (i.e., with A^{r,t) restricted only to the n = band), 
we have the interesting situation of having the continuum gauge-fixing equation 
as well as its discrete realization. For example, in the Coulomb gauge and 
with the restriction to O*'* band, A(r, t) is a simple continuum interpolation 
of its discrete values , each of which is associated with a link £y , where 
£ij connects the z*'* lattice site to its nearest neighboring site j . The discrete 

— * — * 

version of V • ^ = refers to Kirchofi^'s law 

Y,aij = (1.6) 

j 

which is valid at every i , with j summed over all its nearest neighbors. As will 
be shown by (3.37)-(3.39), (4.10) and (4.14)-(4.15) below, in our formulation the 
discrete Kirchoff law implies the validity of the continuum V ■ A — within the 
one-band approximation, and vice versa. Such equivalent formulas can be fur- 
ther generalized. Take V ■ A — 0: in the continuum, this implies that A equals 
the curl of another vector function, say A = V x / , and therefore A becomes 
invariant under the transformation 1^1 + Vx . We are pleasantly surprised 
to find that, within the same one-band approximation in the Coulomb gauge, 
each of these continuum equations has its exact equivalent discrete form (see 
(4.24)-(4.29) below), giving rise to a family of relations between the continuum 
and the discrete. Of course, a change in the gauge has to induce a change in the 
corresponding Bloch functions, as will be illustrated in Section 3 for QED. 

Similar considerations are also applicable to a generalized Coulomb-like 
gauge where, instead of V • A = , we have at all continuum r 

J{f\f\f')-A{f')d^r' = (1.7) 

with f a linear operator independent of ^^(r, t) . While QED is a relatively 
trivial model, the analysis given in Section 3 provides the necessary tools for the 
noncompact lattice formulation of non-Abelian theories. 
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The extension to QCD in the Coulomb or Coulomb-hke gauge is straight- 
forward since the gauge condition, (1.5) or (1.7), is hnear and there are no 
additional supplementary conditions; hence, the group index is external to the 
band decomposition. The details are given in Section 4. Because the noncom- 
pact lattice QCD formulation is the O^'^-band approximation of the continuum 
QCD, it is possible to regard the lattice QCD coupling constant gi as the exact 
renormafeec? coupling constant of continuum QCD, without requiring the lattice 
size £ ^ (i.e., at a finite ), as will be studied in Section 5. We discuss 
how this approach is related to the conventional renormalization procedure in 
continuum QCD. By taking advantage of the asymptotic freedom of the theory, 
we show that, in terms of the perturbation series in the Coulomb gauge, the 
bare coupling qq and the conventional renormalized coupling gn of continuum 
QCD can be readily expressed in terms of the lattice coupling gi^ in a systematic 
way. 

The extension to QCD in the time-axial gauge is more complicated: While 
the gauge condition A^{r,t) = being a linear equation offers no problem 
in the band decomposition, the generalization of the supplementary "Gauss 
Law" condition on the state vector ((3.19) below) becomes nonlinear in QCD, 
which means that in solving this constraint, there would be additional coupling 
between different bands which must be taken into account. The details will not 
be analysed in this paper. 
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2. MATHEMATICAL PRELIMINARIES 



In this section, we consider only the one-space-dimensional case. 



2.1. Lump Functions 

These functions are defined by 

T (^\ — d,k (2 sin 1)"^ j ^ikx/i f^j, ^ ^^^^ 



2tt fc^ \^ik{x-i)/i for m odd 

where m is a positive integer > and x is the continuous coordinate variable. 
(When m = , Lm{x) becomes the 5-function, which does not serve our partic- 
ular purpose.) As we shall see, these functions are nonzero (shaped like a lump) 
only within m lattice cells, each of length i . When m is even and | a; | > ^ £ , 
we may deform the integration path in the complex A;-plane into a semi-circle 
of infinite radius encompassing the upper, or lower, half-plane depending on 
whether x is greater than ^ £, or less than — ^ £; this leads to zero for (2.1). 
Likewise, when m is odd the integral is also zero for |a; — 1 1 > ^ £ . In the next 
section, by using these lump functions we can construct Bloch wave functions 
that satisfy the lattice locality condition (1.2). From the definition, we see that 
for m even 



dx 



-Lm-i{x) + Lm-i{x + tj (2.2) 



and for m odd 



^ = ^m-i(a:) - L^_i(a; - £) . (2.3) 

The lump functions of lower order, m = 1, 2 and 3, are particularly useful 
in later sections. Because, in the application to gauge theories, there will be 
several different kinds of indices, it is convenient to give these functions special 
names. We introduce 

\ T r \ f 1 for < X < £ 

C{x) ^ Li(x) = |q otherwise, 

A(.) ^ L,{x) = f-N<^ (2.4) 

1. U otherwise 

and 

r i[(3-2|f-i|)2-3(l-2|f-i|)2] for |x- IK I 
S{x) = Ls{x) - <^ 1 (3 _ 2 II - i 1)2 for 1 < I X - II < f 

I otherwise . 
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Here, the notation C{x) is chosen for its being a constant within < x < £ , 
A{x) for being hke a triangle within —£< x < £ and S{x) for its nonzero part 
being a second order polynomial in x (within unit intervals of £). 

In general, the nonzero parts of all Lm{x) are positive polynomials of (m — 
-^yh Qp(jej. within each interval. For m > 1 , Lm{x) is continuous at all 
X. It follows from (2.2)-(2.3) that for m > 2 all dLm{x)/dx are continuous 
functions; for m > 3, all d'^ Lm{x)/dx'^ are continuous, etc. When m is even 
> , Lm{x) — Lfn{—x) ; it peaks at x = , decreases monotonically to zero at 
\x \ = ^ ^ cind remains zero for \ x \ > ^ £ . When m is odd > 1 , Lm{x — |) = 
Lm{—x + f ) ; its maximum is at x = | , decreases monotonically to zero at 
\x — ^ \ = Y ^ ^^'^ remains zero for | a; — | | > ^ ^ • 

Another useful property is, within the range < x < £ , 

m/2 

Lm{x — j£) = constant for m even (2.5) 

j=-(m-2)/2 

and 

(m-l)/2 

Lm{x — j£) — constant for m odd, (2-6) 

j = -(m-l)/2 

which can be readily proven by noting that their derivatives are both zero, on 
account of (2.2)- (2.3). 



2.2. Block Functions 

Given an m , the corresponding (unnormalized) O*'^ band Bloch function 
is given by 

oo 

T{K\x) = e'^'L^{x-j£) (2.7) 



with 



which satisfies 



Because 



e = K£, (2.8) 

16*1 < TT. (2.9) 

J'{K\x + £) = e'^T{K\x), (2.10) 
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the passage from an infinite lattice to a finite one of size Ni and with the 
periodic boundary condition can be readily derived by requiring 



KNl = 27r X integer. (2.11) 

This convention will be used subsequently. It follows from (2.5) and (2.6) that, at 
X = , !F{K \x) is a constant at all x . Furthermore, the sum A^~^ YIk e~^^^^ 
xT[K I x) = Lm{x — j£) satisfies the lattice locality condition (1.2). 

To construct the higher-band Bloch functions, we follow one of the simple 
procedures given in I. Let 

iPpix) = {N£)-^ e^P^ (2.12) 
be the usual plane- wave functions with pN£ = 2tt x any integer. Resolve 

P = if + ?f (2.13) 

where n = •• , -2, -1,0, 1,2, •• . At a K given by (2.8)-(2.9), write (2.12) as 

^^{x) = (Ar£)-^e^^^e*2W^. (2.14) 

The complete set of orthonormal Bloch functions, given by {fn{K \ x)} , can be 
obtained by requiring, for n = 

MK\x) = CmHK\x), (2.15) 
with Cm the normalization constant, and for n ^ 

fn{K\x) = Ipnix) — an (2.16) 

1 + ao 

where (for all n , including 0) 

an = / fo{K\xy iPn{x)dx (2.17) 
Jo 

and * the complex conjugation. It follows then 

/ U{K I xYU'iK' \x)dx = 5kK' Snn' (2.18) 
Jo 



where 5kk' and 5nn' are Kronecker symbols. By using (2.1) and (2.7) we see 
that 









1 Cm 1 



2 2"^(1 -cos6')" 
{9 + 27rn)2^ 



(2.19) 



Since X]n=-oo I | = 1 , we may write (2.19) as 



|onP = (6' + 27rn)-2"^ 
Thus, when ^ — > , 



-2m 



.n — — oo 



-1 



ao 



l_2^2my^27rn)-2- = I - 

' (2m)! 

n=l ^ ' 



and, for n , I a„ p ^ 6^'^{2TTn) ^"^ where 5^ denotes the BernouUi 
numbers with Si = | , ^2 = ^ , = , etc. 

We note that, given J^{K \ x) , there are many equally simple ways to con- 
struct the complete set of Bloch functions, as discussed in I. 



Remarks: 

1. When m = 0, L^(x) becomes 6{x/i) , as mentioned before; therefore, the 
expansion (2.7) reduces to the usual expression for a lattice of discrete points. 

2. In (2.19), we see that for nonzero n and m the factor {9 + 27rn)~^'^ is 
always less than [ (2 | n | — 1) tt]"^™ . In D-dimension, (2.19) will be changed 
into a product of D such factors, one for each space dimension. Thus, the 
higher band Bloch functions are very close to the corresponding Fourier series 
of high wave numbers (see, e.g., Eq. (3.67) below). 



2.3. Scalar Field (analysed in terms of Lm{x) for m = 2) 

As an example, consider the Lagrangian for the free scalar field (f){x, t) , in 
units fi — c— 1: 



Nl 



dx dx 



dx 



(2.20) 



where fj, is the mass, ^ denotes the Hermitian conjugation and the dot indicates 
the time derivative. The restriction to the zeroth band n = by using the lump 
function of order m = 2 is identical to the expansion 



(2.21) 



where j goes over the N lattice sites and A (a;) is the triangular function given 
by (2.4). At the j^^ site x = j£ , = 0^ and the expansion (2.21) gives a simple 
linear interpolation, converting the set of discrete values { (pj } to a continuous 
function. (The use of piecewise flat but continuous functions such as A(a; — ji) 
for lattice formulation is the basis of discrete mechanics.'*) 

We may substitute (2.21) into JC and derive the corresponding lattice equa- 
tion by setting to zero the variational derivative of the action J Cdt with respect 
to (t) . The resulting solution is 



K 



Ne{2 + cos0) 



(2.22) 



where Q%{t) is the normal mode satisfying 



with 



6(1 - cos^) 



£2(2 + COS ^) 



(2.23) 



(2.24) 



and 9 given by (2.8)-(2.9). Correspondingly, the lattice approximation (2.21) 
becomes 

cf>{x,t) = J2QK{t)MK\x) (2.25) 



K 



where the normalized /o (K \ x) is 



fQ{K\x) 



N£{2 + cose) 



(2.26) 



Clearly, the full expansion in terms of all bands 



0(x,t) = Y,Y.QK{t)fn{K\x) 



n K 



(2.27) 
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yields the exact continuum result. The details are given in I. 
Remarks. 

1. For K within the Brillouin zone (i.e., |^| < tt), the above frequency 
oj{K) = \J'w'^{9) + iJ? is always higher than the exact continuum result 
E{K) — K"^ + fj? . This is expected since the expansion (2.21) may be re- 
garded as a constraint imposed on the continuum problem; the corresponding 
optimal spectrum can only serve as an upper bound. 

2. Instead of A(x — j£) , in the expansion (2.21) for (p{x,t) we may use another 
lump function Lm{x — j£) , but with m > 2 , as a constraint on the continuum 
field. (See, e.g., (2.28) below.) However, this does not work if m is < 2 : When 
m = 1 , Li{x) — C{x) and its derivative is d{x) — d{x — £) , in accordance with 
(2.4); therefore, the application of the corresponding 0*'*-band Bloch function 
to the continuum scalar field Lagrangian C would lead to oo . (However, the 
m = 1 lump function can be used as a link function for the gauge field, as will 
be discussed in Section 3.) When m = 0, Lo{x) = 5{x/£) ; the substitution of 
the corresponding expansion of (f){x) into the continuum Lagrangian C would 
clearly also give oo . 

3. In this context, it may be useful to examine the usual discrete lattice of 
N separate points. The corresponding Lagrangian Cd for a discrete field ipj 
would be, instead of (2.20), 

N 

In terms of its normal modes we have, in place of (2.22)-(2.24), 

K 

qk = -uJd{KfqK 

with 

ujd{K) = [2(1 -cos^)£-2 + //2]^ . 

Note that the frequency uJd{K) is /ower than the continuum result E{K) , which 
underscores the difference (for ^ 7^ ) between a discrete lattice Lagrangian Cd 
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and its continuum analog L. When £ — > , '^{K) and LOdiK) both approach 
E{K) with the deviations, uj{Kf - E{Kf and uJd{Kf - E{Kf , 0(^^) . 

However, when £ 7^ , the usual discrete lattice Lagrangian is more an 
abstract model of the continuum theory, rather than a systematic approximation 
of it. This is different from the direction that we follow, and the difference 
becomes especially prominent when we come to QCD (Sections 4 and 5). In 
Appendix A, we give a detailed comparison between these two approaches, which 
may help to explain the underlying reasons for constructing these special lump 
functions and for the particular way that we formulate our lattice approximation. 



2.4.- Link Functions 

For m odd, Lm{x — ji) peaks at x = {j + , which is the mid-point of 
the link connecting the j*'^ and (j + 1)*'^ sites, located at x = ji and (j + 1)£ . 
It is convenient to associate these functions with links, as will be discussed in 
later sections; therefore, the lump functions of odd order will also be called "link 
functions" . Here, in one space-dimension, this difference is not an essential one 
since the mid-point of any link is also the lattice site in the dual structure. We 
shall illustrate the use of these link-functions by considering the special case 
m = 3 . 

For the scalar field (j){x, t) , instead of (2.21), write 

<t>{x,t) = J2xj{t)Six-je) (2.28) 

where S{x) = L-:i{x) is given by (2.4). Substituting this form into the La- 
grangian (2.20) and varying Xj{^) 1 we find, similar to (2.22)-(2.26), 



K 



30 



2 



go^^^^ ^2.29) 



7V£(16+ 13cos^ + cos2 e) 

where Q%{t) is the new normal mode satisfying 

Q% = -iv\e) + fi^)Q% (2.30) 



with 



2,„, 20(1 -cos ^) (2 + cos ^) 



^ - £2(i6+i3cos^ + cos2^) ^^-^^^ 
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and 6 given by (2.8)- (2.9), as before. Correspondingly, (2.28) becomes 



^{x,t) = Y,QK{t)9o{K\x) (2.32) 



K 

where 



gQ{K\x) 



30 



iV£(16 + 13cos6' + cos2 9) 



e'^^+'^^^ S{x - jt) . (2.33) 



For notational clarity, we denote the normalized zeroth-band Bloch function for 
the link functions by gQ{K \ x) , instead of fo{K \ x) . 

From (2.4), we see that 5(|) = | and S'(-|) = S{f) = i. Thus, the 
value (pj{t) = (p{x,t) at the mid-point of each link (i.e., at x = {j + ^) £) is 
related to Xj(t) introduced in (2.28) by 



(6x,+Xi-i+Xi+i). (2.34) 



The inverse relation is 

N 

Xj = J2Y1 4[iV(3 + cosa)]-' 0,., (2.35) 

a j' = l 

where 

a = 2tt n/N 

and its sum extends over n= 1,2, ■■ ,N . Both Xj 0^ satisfy the periodic 
boundary condition Xj+N = Xj ^j+N — 4'j • 

In the usual lattice calculation, one tends to use (f)j as the primary variables. 
To see how closely correlated the corresponding Xj ^^^i us consider the 
special case N ^ oo and 



= 1 1 ^orj = (2.36) 
1 otherwise 



^ 1 otherwise 
The corresponding Xj is 

Xj = V2 Zl ^' I ^ {-y V2 e-^-^^^ I ^ I (2.37) 

where Z = —3 + 2-\/2 . Thus, Xj decreases rapidly with increasing 

13 



By following the steps (2.16)-(2.18), we can start with the 0*'* band link- 
function go{K I x) and construct a complete set of orthonormal Bloch functions 
{gn{K\x)}, satisfying 

/ gn{K\x)* gn'{K' \x)dx = 6kK' Snn' ■ (2.38) 
Jo 

The expansion of (f){x, t) in terms of all bands would, of course, give the correct 
continuum solution. 

For 1^1 < TT , like w'^{9) , v^{d) is larger than the exact continuum eigen- 
value K"^ = {Q jiY • This is again expected since the expansion (2.28) may also 
be regarded as a constraint imposed on the continuum problem. In addition, we 
have 



2(1 -COS 6*)^ (8 + 7 COS 6*) 
a relation that will be of use later on. 



2.5. Generalized Lump Functions 

The generalized lump functions are defined by 

L (x) ^ r - (2 ^"4)"' ('^•"4)" ^ r e^^-/^ for m + n even 

""'"^ ^ i-oo27r \^ikix-^)/£ form + nodd 

(2.40) 

where m and n are both positive integers. When n = 0, Lm,n{x) reduces 
to Lm{x) given by (2.1). These generalized functions are nonzero only within 
m + n lattice cells. Similar to (2.2)-(2.3) and (2.5)-(2.6), we have for m + n 
even 



^ ^ -^m,n(a;) = -Lm-l,n{x) + Lm-l,n{x + £) , (2.41) 
(m+n)/2 

Lm,n{x — j£) = constant 

j=-(m+n-2)/2 



and for m + n odd 
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(■ ^ L^,n{x) = Lm-l,n{x) - Lm-l,n{x - i) , (2.42) 

(m+n-l)/2 

Lm,n{x — j£) = constant. (2.43) 

j=-(rn+n-l)/2 

The following example for m + n = 4 illustrates the variety of ways that 
we may interpolate from a discrete set of lattice values {4>j} to a continuous 
function (j){x) in the continuum. Introduce a lump function of mixed orders: 



L{x) = 3L4,o(x)-2L3,i(a;). (2.44) 



We find, 



|2 



+ for \x\<e 

L{x) = <j i(i_M)(2-M)2 iove<\x\<2e (2-45) 

for I a; I > 2£ . 

Because L(0) = 1 and L{±£) = , the expansion 

0(x) = <j>^ L{x - jl) (2.46) 

j 

gives (l){xj) = at a/Z lattice sites 0;^ = , even though each (f)j influences the 
value of 0(a;) over four nearby lattice cells (instead of two, as in the expansion 
(2.21) using A(a;) ). The above lump function L{x) is not always positive; it is 
1 at X = , decreases to at x = ±£ , reaches a minimum -2/27 at a; = ±|£, 
then increases to at x = ±2£ and remains for | a; | > 2£ . The corresponding 
normalized 0*'* band Bloch function is, for < x < £ , 



k{K\x) 



70 



X 



_N£{<o2 + 17 cos 6* - lOcos^ 6 + cos^ 6* J (2.47) 
[e-'^ L{x + £) + L{x) + e'^ L{x - £) + e^'^ L{x - 2£) ] 



with Q given by (2.8)-(2.9). From (2.47), one can construct a complete set of 
orthonormal Bloch functions { /„(i^ | a;) } , as before. 
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3. QED (Free Field) 



The Lagrangian density of a free electromagnetic field {f, t) is 




(3.1) 



where F^jy = — -g^ , and the repeated subscripts jj, and z/ are summed 
over 1, 2, 3 and 4 (i.e., x^y.z and it). We follow the standard procedure of 
deriving the continuum quantum Hamiltonian by first fixing the gauge. The next 
step is to choose an appropriate complete set of orthonormal c number Bloch 
functions, in terms of which the quantized field operators can be expanded. The 
restriction to the 0*'* band Bloch functions constitutes the lattice theory in the 
present approach, as in (2.27) and (2.25). A critical point for the gauge theory 
lies in being able to preserve the continuum gauge-fixing condition (valid at all 
continuous space-time coordinates) within the lattice approximation. Because 
of (1.2), the lattice gauge theory maintains the local character of the original 
field theory. The compatibility between the continuum gauge-fixing and the 
band decomposition makes it possible to carry out the transition from lattice 
to continuum in a straightforward and systematic way. As we shall discuss, the 
result is a new lattice gauge formulation, which retains the noncompact nature 
of the continuum gauge theory, and thereby may serve as its first approximation. 

3.1. Time- Axial Gauge 

In the time-axial gauge. 



Embed a cubic lattice oi M = sites in the continuum. Consider first the 
unit cubic cell (of volume ) which has the origin (0, 0, 0) and the lattice-site 
(1,1, 1)£ as two of its eight vertices. Altogether the cell borders on 12 links. 
Assign the values of Vx on the four x links (all parallel to the x-axis) to be the 
constants a, a', a" and a'" . Likewise, on each of the y and z links Vy and Vg 
assume constant values 6, •• , b'" and c, •• , c'" , as shown in Figure 1. Inside the 
cell, let be a multilinear function of y and z : 



V4 = iVo = 0. 



(3.2) 



Vx = a+ {a' - a) ^ + {a" - a)^ + {a'" - a" - a' + a) 



yz 



(3.3) 



£2 ; 



16 



similarly, 



Vy = b+{b'-b)^^+{b"-b)^^+{b"'-b"-b' + b)'^ (3.4) 



and 



V. = c+{c'-c)^^+{c"-c)^^ + {c'"-c"-c' + c)^. (3.5) 
The magnetic field B = V xV inside the cell is given by 



= i 



-1 



(l-j) ic" -c-b' + b) + ^^ {c'" -c'- b'" + b") ] , (3.6) 



and similar expressions for By and Bg . Thus, on the (y, z) plaquette at 
X = 0, B^ is a constant. From Figure 1, we see that this constant value 
£~^{c" — b' — c + b) is equal to times the algebraic sum of link- values 
along the counter-clockwise direction, with the signs determined by the arrow 
directions shown in the figure. Likewise, on the (y, z) plaquette at x = £ , B^ 
is again a constant, i~^{c"' — b'" — c' + b") . Inside the cell, Bx depends linearly 
on X . Identical considerations can be applied to By and B^ . 

Next, we consider other cubic cells, each of which (say, the j*'* cell) can 
be arrived at through a space-translation, moving the corner of the above unit 
cell from the origin (0,0,0) to the j^'^ lattice site at fj = (ji,j2,i3)^- At the 
same time, replace the link- values a,b and c by aj,bj and Cj . The entire 
Vp-function in the continuum space can then be obtained by following the same 
multilinear interpolation procedure (3.3)-(3.5). This will lead to the 0*'* band 
approximation in the new lattice gauge formulation. In this way, the aj , bj and 
Cj are defined once and for all on each link, so that in passing from one cell to 

— * 

another, there is no discontinuity of either the tangential components of V or 
the the normal components of B . (See (3.6) above.) 

As we shall see, by using the lump and link functions introduced in the 
previous section, we can cast the result in terms of Bloch functions consisting 
of products of (2.7). We write 

Va = J2 (la{K)g{Ka\Xa) J] f{Kb\xb) (3.7) 

where the subscript a or 6 denotes the space-components 1,2 and 3 (or x^y 
and z), K denotes the Bloch wave number, which is to be summed over the 
Brillouin zone with its components Ka satisfying (2.8), (2.9) and (2.11); i.e. 
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Ka N£ = 27r X integer (3.8) 

— * — * 

and the components of ^ = Ki bounded by 

-TT < da < TT . (3.9) 

The functions g{Ka \ Xa) and /(i^b I ^h) ai'e 

g{Ka\xa) = E e^"'" ^(^a - n£) (3.10) 

and 

^ ^ ' ' n= — OO 

where C(a;) and A{x) are the lump functions Li{x) and L2{x) introduced in 
the previous section. By referring to (2.4), we see that (3.7) reproduces (3.3)- 
(3.5). 

Consider again the j*^ cubic ceU. By using (2.4), (3.7) and (3.10)-(3.11) 
we can write down immediately the values of Vx,Vy and Vz on the three links 
connected to the lattice site at fj . Thus, 



= E {j^y (jtLs;)' -(^)- 

K 

1 1 

bj - y ^ e'^-''' ( ^ ^ ' ( ^ ^ ' q2{K) , (3.12) 

^ 4^ V2 + COS^3y \2 + 00801 J ^ ^ ^ ' 

K 

where Q, = {N£)^ is the total volume. (In this connection, it may be helpful 
to compare the conventional discretized Abelian lattice gauge theory to the 
present formulation based on using Boch wave functions, of which the zeroth- 
band ones are simple continuum interpolations of discrete values aj , bj and cj . 
A discussion is given in Appendix A.) 

Substituting (3.3) and (3.7) into (3.1), we find the Lagrangian (in the 0*'* 
band approximation) to be 



18 



Cd^r = [qHK)q{K)-qHK){w^ -ww)q{K) 

K 



(3.13) 



where q{K) and w = w{$) are 3x1 column matrices whose components are 

QaiK) and 

with 9 = Ki as before. Sometimes these column matrices will also be repre- 
sented by vectors q{K) and w{9) . In (3.13), denotes the transpose and 



^ = |:-(«T = E^^, (3.15) 



which is the three-dimensional generalization of (2.24). Throughout the paper, 

I — * 

^ ^ denotes the sum over half of the JC-space within the Brillouin zone. The 
corresponding Hamiltonian is 

H = J2 [pHK)p{K)+q\K){w^ -ww)q{K) (3.16) 

K 

with p{K) = q{K) satisfying the equal-time commutation relation 

Both the Lagrangian and the Hamiltonian are invariant under the time- 
independent gauge transformation 

aj ttj + (pj' - (pj , 

bj — > bj + (f)j" — , (3-17) 

Cj Cj + C[)yn - (Pj , 

where the subscripts refer to the sites 

j = (ii,i2,i3), / = (ji + I,i2,i3) , 
j" = ih, 32 + 1,33) and = {ji,j2,h + l). (3.18) 
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This property can be readily established by noting that V is unchanged since 
the transformation is independent of time and B = V xV is also not affected, 
on account of (3.6). As in the standard treatment in the time-axial gauge, the 
above symmetry necessitates the following constraint on the state vector | > : 



Jj\>= (3.19) 

where Jj is the generator associated with the transformation (3.17). There 
are altogether J\f — 1 such constraints. (Note that the special case 0i = 02 = 
• • • = produces no change in a-,-, hj and Cj ; therefore (3.17) for all j only 
imposes M —1 constraints.) 

In the Hilbert space representation in which the basis vectors are chosen to 
be the eigenstates | > of the operators qa{K) , the state vector takes the form 
< q\ > and (3.19) becomes 

Jj<q\ > = (3.20) 
with Jj given by the differential operator 



3 



Hence, (3.20) is equivalent to, for K ^0 , 

3 Q 

E^«^7^ <q\>= 0- (3.22) 



At a given K , introduce a right-hand basis of three ortho normal vectors 
ei{K) , e2{K) and w{K) = w/\/w^ . Define the generalized coordinate for the 
longitudinal mode to be 

Qe{K) = w{K) ■ q{K) (3.23) 
and that for the transverse mode 

Qt{K) = et{K)-q{K) (3.24) 

where t = 1, 2 and q, w are the vectorial forms of the column matrices q and 
w. The constraint (3.22) implies < q\ > to be independent of all Qi{K) . 
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Thus, in the 0*'* band approximation, the spectrum of the free electromagnetic 
field is given by two transverse modes for each K in the Brillouin zone with 
w'^{9) as the square of the energy. (Use of the subscript I notation to denote 
"longitudinal" occurs only in this subsection and Appendix A.) 

It is straightforward to construct higher band Bloch functions through 
(2. 16)- (2. 18). Expanding in terms of all bands, the spectrum changes from 
w'^{9) to the exact continuum result = Q'^ jP' for K within the Brillouin 
zone, as expected. 

Remarks: 

The inverse relation of (3.23)-(3.24), 

q{K) = w{K)Qi{K) + J2^t{K)Qt{K), 

t 

enables us to decompose (3.7) into a sum of two parts: 

V{f,t) = Vt{f,t) + V'{f,t) 
with their components given by 

iVi)a = X;^(^)aQK^)^(i^a|a;a) J] fi^blxb) 

and 

From (2.2) and (2.4), we have 

^ /\{x) = -C{x) + C{x + tj (3.28) 
ax 

which leads to the following relation between g{K\x) and f{K\x)^ defined by 
(3.10)-(3.11): 

I ^ f{Ka \Xa) = i Wa g{Ka \ Xa) (3.29) 

with Wa given by (3.14). Thus, V x = at all points in the continuum 
space, since 
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(3.25) 
(3.26) 

(3.27) 



Ve{f,t) = Vxir.t) (3.30) 

where 

X = ^Xk Y{f{Ka\xa) 

K « 

and (3.31) 

= -iQi{K)/V^. 

On the other hand, because dg {Ka \ Xa) / dxa is very different from f{Ka \ Xg) , 
one can readily verify 

V ■ V' ^ 0. (3.32) 

Indeed, should we require V' to be divergence free, the roles of / and g would 
have to be interchanged in (3.29), as will be discussed in the next section. 

The state vector < g | > introduced in (3.20) is a functional of V (i.e., of 
and V') . The constraint (3.22) makes < q\ > to be a functional only of 
V'{f, t) , independent of the longitudinal component Vi{r, t) . However, V-V ' ^ 
reminds us that the familiar difference between the time-axial gauge and the 
Coulomb gauge persists throughout the band structure. 

We note that in the decomposition (3.25), the component et{K) Qt{K) 
is, of course, perpendicular (i.e., transverse) to the longitudinal component 
w{K)Q(,{K) . Because K lies within the Brillouin zone, this is only the lat- 
tice transversality condition; in contrast, (3.30) and (3.32) refer to relations in 
which r is a continuum coordinate. 



3.2. Coulomb Gauge 

Let Afj_{f, t) be the electromagnetic field in the Coulomb gauge. Its spatial 
part A{f.,t) is divergence-free; i.e., 

V ■ A = 0. (3.33) 

The free-field Hamiltonian is 

H = I ^{n^ + B^)d^r (3.34) 
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where S = V x yl denotes the magnetic field and — H is the electric field which 
also satisfies the transversality condition 

V • n = 0. (3.35) 
In addition, we have the equal-time commutation relation 



(3.36) 



The central problem that confronts us is the construction of the appropriate 0*'* 
band Bloch functions, so that (3.33) and (3.35) can be valid at all continuum 
points in the lattice approximation. 

As mentioned at the end of the preceding section, in order to have a 
divergence-free field, the roles of / and g in (3.29) have to be interchanged. 
Thus, we replace the link function (3.10) by (2.33). We write, for the O*'^ band 
function expansion. 



K 



and 



^a{f,t) = J2pa{-K)Fa{K\f) 



K 



where 



Fa{K\r) = Cs{ea)g{Ka\Xa) J] ^2 (^b) ^(i^6 | ^^b) 

in which, similar to (2.7), but with m — 2 and 3, 



(3.37) 



(3.38) 



and 



(3.39) 



g{Ka\xa) = e''^'- J2 e'^'-S{xa-ne). 



The factors 



C2{0a) 



Ne(2 + COS Oa) 
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and 



(3.40) 



cM = 



30 

Ne{16 + 13 cos Oa + COs2 Oa) 



refer to the normalization constants Cm in (2.15) with the subscripts m = 2 
and 3 denoting the order of the lump functions in and Q . We note that apart 
from normalization factors, and Q are the same functions as / of (3.11) and 
5^0 of (2.33). 

From (2.3)-(2.4), it follows that 



dS{x) 
dx 



= A{x)-A{x-£) 



and therefore 



d 9 

— g{Ka I Xa) = 2i sin y J^{Ka I Xa) 



The transversality conditions (3.33) and (3.35) become simply 



and 



q{K) ■ v{K) = 



p{K) ■ v{K) = 



where the components of v{K) are 



Va[K) = ;^ smy 



10 (2 + cos Oa) 



Because 



we have 



16+13 cos 9a + cos^ 9a 
Fa{K\rY = Fa{-K\r), 

,{K) = Qai-K^ and p«(K) = Pai-K)^ 



(3.41) 
(3.42) 



(3.43) 



(3.44) 

(3.45) 
(3.46) 



— * 

At a given K , form a right-hand basis of three real orthonormal vectors 
ei(^), e2iK) and v(K) = v{K) / ^ v{KY so that 



ei{K) = e2{K) X v{K) 



(3.47) 
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Thus, if e2{-K) = e2{K) then ei{-K) = -€i{K) , since v{-K) = -v{K) , in 
accordance with (3.44). Write 



and 



(3.48) 



,XK) = 5^et(i?)„Pt(i?) 



where the sum is over t = 1, 2 . The commutation relation (3.36), when applied 
to the O*'^ band components of 11 and A , gives 



Pt{K),Qt,{K') 
or in its equivalent form 

p„(i?),g„,(K')] = -i5j^j^,(5aa' -Va{K)Va'{K)y 



(3.49) 



(3.50) 



In the lattice approximation (3.37) and in terms of the column matrices 

(Pi{K)\ (qi{K)\ 
p{K) = p,{K) and qiK) = q,{K) , (3.51) 

\P3{K)J \q3iK)J 

the Hamiltonian (3.34) becomes, on account of (3.28) and (3.43), 



H = ^[p^p + q^ Aq' 



(3.52) 



K 



where the matrix A = (^Aab{K)j is diagonal, with its diagonal elements 



given by (3.15) and, similar to (2.39), ba = w\ — \ i.e.. 



£2 (2 + COS 6'a) (16 + 13 COS Oa + COs2 ^„) ' 



(3.53) 



(3.54) 



— * 

As before, Oa = Ka £ . At a given , although A is a 3 x 3 matrix, because of 
the orthogonality condition (3.43), there are only two transverse modes whose 
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eigen-frequencies a;^ = o;^ are determined by the secular equation (cf. Eq. 
(2.16) of I) 



+ 



All - A22 - A33 - a;^ 



0, 



(3.55) 



The solution is 



=w'-^ [vl{52 + h) + vl{6:, + 5i) + ^2(^1 + ^2) 



[ 8l{vl + vlf + 5l{vl + vlf + 5l{vl + vlf 



2^2 



(3.56) 



where 



V r^l2 ^ 20(l-cos^J(2 + cos^J 

Z^^«U 2^ £2ne^i3cos6'a + COs2 6'a) 

a=l ^ ' 



(3.57) 



0=1 



Because 5^ > , each of the diagonal elements Aa^ lies within v'^ and . 
The two eigenvalues have to be sandwiched between An , A22 and A33 ; 
consequently, they must also satisfy 



2^2^ 2 
V < UJ_^ < w . 



(3.58) 



Using the spherical coordinates v\ = v sin a cos /? , V2 = v sin a sin /? and V3 
V cos CK , we may express (3.56) in an alternative form: 



1 



All + A22 ± V (^11 - ^22)2 + 4A 



12 



(3.59) 



where 



A22 — 



and 



All = —Si sin^ /? — ^2 cos^ /? , 
-{5i cos^ P + 82 sin^ /?) cos^ a — 6^ sin^ a , 

A12 = {Si — 82) sin/5 cos/? cos a. 



In the long wavelength region | 9a \ « 1 , deviates from w'^ of (3.15) by 
0{9'^) , and (3.59) becomes 



a;. = w 



24^2 



S2 



3S3 
Si 



■S3 

S2 

Si 



S2 



9S3 
Si 



+ 0{e^) (3.60) 
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where si = 6^ + + , S2 = + ^| + $1 dl and S3 = Ol Ol ^3 • 

The higher-band Bloch functions can be constructed through a generaliza- 
tion of (2.16)-(2.18). (See (3.73)-(3.79) below). As we shall see, these higher- 
band Bloch functions are very close to the corresponding Fourier series of high 

— * 

wave numbers. By expanding A in terms of all bands, one changes the spectrum 
from the lattice approximation uj\. to the continuum result. 

Remarks: 

1. By using (3.48), the 0*'* band expansion (3.37) of A in the Coulomb gauge 
can be written as (with the time variable suppressed) 

2 

it *=i 

where F^{K | r ) is a vector function whose components are 

F\K\r)a = et{K)aFa{K\f)- (3-62) 
Because of (3.38)- (3.42) and (3.47), for any K within the Brillouin zone 

V ■ F\K\f) = (3.63) 
everywhere in the coordinate space. Thus, we have 

F\K\f) = WxG\K\f) (3.64) 
with the components of G^{K \ f) given by 

2 

t' = l 67^0 

(3.65) 

and {rf^') denoting the 2x2 antisymmetric matrix whose elements are r]^^ = 
rf''^ = and ■q^'^ = —rf^ = 1 . Although the independent variables in the 
lattice approximation are restricted to the set {Qt{K) } , the 0*'* band Bloch 
functions F*{K \ r) satisfy the differential relations (3.63)-(3.64) at all r in the 
continuum. 
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2. From (2.1) and (2.4), we see that the Fourier series expansion of (3.38) is 

Fa{K I f ) = e'^-^ ^ Ca(^| m) e^'^^^'^l^ (3.66) 

where 

TTT-a (the components of m) are all integers, and C2, C3 are given by (3.40). 

3. To construct the higher-band Bloch functions in the Coulomb gauge, we 
generalize the steps (2. 12)- (2. 16): Define 

i)^\r) = n-U\p)e'P-^ (3.68) 

where t = 1,2 denote two sets of divergence- free functions, with and 
both real unit vectors perpendicular to p and to each other. Resolve 

p = K +{2tt fi/t) (3.69) 

with the components of n being any integer. At a given K within the Brillouin 
zone, write (3.68) as 

^-*(f) = n-5e*(n)e^^-^^+*2'^^-^/^ (3.70) 
Introduce the 2x2 matrix Tfi whose elements are 

T'J'{K) = J F\K\r)-ifn''{r)*d^r (3.71) 

where F*{K \ r) is given by (3.62). The complete set of orthonormal divergence- 
free Bloch functions {Fi{K\r)} can be constructed by requiring, for n = , 

F^{K\r) = F\K\r) (3.72) 

and for n 7^ 

2 

Fi{K\f) = ^n\r)-Y, [(M^)"'F*'(i?|f) + (iV^)"Vo*'(r)] . (3.73) 

t'=i 
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— * , — *./ 

The orthogonality between Fq and Fi yields the 2x2 matrix relation 

Tn-Ml-ToN^ = 0, (3.74) 

which together with the orthonormality between Ft and Ft for nonzero n 
and m give 

MnM^.^NnN^. (3.75) 

The solution of (3.75) is 

Nn ^ Mnu (3.76) 
with u an arbitrary 2x2 unitary matrix satisfying 

uu^ = 1. (3.77) 

By using (3.74), we have 

and (3.78) 

Nn = 4{1 + uTS)-'u. 

It follows then (for all n and rh , including 0) that the orthonormality relation 

/ d\Fi{K\r)* ■ Fi{K'\T) = SnrnStt'Sj^g,; (3.79) 

in addition, every member of the set {Ft} satisfies the divergence- free condition 

V-Fi{K\f) = 0. (3.80) 
If one wishes, the arbitrary u can be set to 1; (3.78) becomes simply 

= = T^(l + 4)-^ (3.81) 

Because of (3.67), the matrix elements 

3 

= ^e*(i?)a-e*'(n)„C„(^"'|n) (3.82) 
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are all real. The expansion of A{f,t) and n(r, t) in terms of all bands, 

A{r,t) = J2 Qn{K)Fi{K\r) 

n,K,t 

and (3.83) 

n,K,t 



changes the O*'* band spectrum from the lattice approximation uj^.{9) of (3.56) 
to the exact continuum result = {0 /£)'^ within the Brillouin zone and, for 

1 2 



fi^O, K + {27Tn/£) 



outside the zone. 



4. The gauge- fixing conditions, Vo{f,t) = in the time-axial gauge and 

— * — * 

V- A(f, t) = in the Coulomb gauge, are valid at all f in the lattice approxima- 
tion. This necessitates the difference between the 0*'* band functions in these 
two gauges. Of course, in either gauge, the Bloch functions of all bands are 
complete. The transformation matrix between these different band structures 
can therefore be readily derived. 

5. Similar considerations are also applicable to a generalized Coulomb- like gauge 
where, instead of V ■ A — Va Aa — , we have^ at all continuum r 

j{r\Va\f') Aa{f')d'r' = 0, (3.84) 

with Ta a linear operator independent of A . Accordingly, we have to create 
a new class of lump functions of different orders m , for which relations like 
(2.2)-(2.3) remain valid, but with the differential operator d / dxa replaced by 
Ta ■ Details will be given elsewhere. 

The extension to QCD in a Coulomb or Coulomb-like gauge is straightfor- 
ward since the group index is external to the band decomposition, as will be 
discussed in the next section. However, as mentioned in the Introduction, the 
extension to QCD in the time-axial gauge is more complicated: While the gauge 
condition ^4(r, t) = being a linear equation offers no problem in the band 
decomposition, the generalization of the constraint equation (3.19) 

Jn \ > = 
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becomes nonlinear in QCD. That means in solving this constraint, there would be 
additional coupling between different bands which must be taken into account. 

In the next two sections, we will discuss QCD, but only in Coulomb gauge. 
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4. QUANTUM CHROMODYNAMICS 



Let ^^(r, t) be the field in tlie Coulomb gauge, with the subscript a = 
1,2,3 denoting its spatial components and the superscript i = 1,2, ■■, — 
1 its SU{n) index. The conjugate momentum is V^{r,t) . Both satisfy the 
divergence-free condition 



V ■ A\f,t) = V • V\f,t) = 



and the equal-time commutation relation 



(4.1) 



Vi{f,t), Ai,{f',t)\ = -,5^^(5,,-V-2V„Vb)5^(f-f') 
The continuum QCD Hamiltonian (without quarks) is^ 

'^ = \ j [J-'rijvi + BiBi]d'r + nconi 

where the Coulomb interaction is 



(4.2) 



(4.3) 



with the charge density 

the color magnetic field Bf^ is given by 



Va denotes the covariant derivative matrix whose elements are 



ya - go J 



J is the the Jacobian^, or the Faddeev- Popov determinant^ 



J = det I Va^al, 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 
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Sabc is the usual antisymmetric tensor in the three-dimensional space, J*"^" is 
the antisymmetric structure constant of the SU (n) group algebra and is the 
bare coupling constant, related to the renormalized coupling g by 



go = g + Sg, (4.9) 

as will be discussed later. All the above repeated indices are summed over (but 
not for the equations below). 



4.1. 0^^ Band Expansion 

The set Fi{K\f) , given by (3.72)-(3.73), forms a complete orthonormal 
functional basis that satisfies {V ■ Fi{K\f)} = . Consequently, we can al- 
ways expand the field operators A^{f,t) and V^if^t) in terms of these Bloch 
functions, similar to (3.83). In the O*'^ band approximation, we restrict the 
expansion only to n = 0. For clarity, the O*'* band approximations of A^{f^t) 
and V^{f,t) will be denoted by A^{r^t) and n^(r, t) , whose components are 
given by, analogous to (3.37), 

Km = Y.^a{K)F^{K\f) (4.10) 

K 

and 

nl(r,0 = Y^pi{-K)Fa{K\f) (4.11) 

K 

where Fa{K\ r) is given by (3.38). As in (3.43) and (3.50), the divergence-free 
condition (4.1) requires 

3 3 

= E/aW^a(i?) = 0, (4.12) 
a=l a=l 

and the commutation relation (4.2) leads to 

pi{K),qi,{K') ] = -i5'''5^^, (5aa'-Va{K)Va'{K)) . (4.13) 



The substitution of (4.10)-(4.11) into (4.3)-(4.4) gives the noncompact lattice 
formulation of QCD; in addition, g^ should be replaced by g^ , the lattice cou- 
pling constant, as will be discussed in Section 5. (See (5.1)-(5.6) below.) 
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Because of (3.38)-(3.39), the Bloch function Fa{K \ r) satisfies the lattice 
locality condition (1.2). To appreciate its consequences, we express (4.10), the 
expansion of A^{f,t) , in an equivalent form: 

Ai{r,t) = J2^'j(t)S{x-h£)A{y-h£)A{z-h£), 

j 

K(^^ t) = 11 ^(^ - h^) S{y - j2i) Aiz - 3d) (4.14) 



and 



4(f,t) = ^c^(t)A(x-ji£) A(y-j2^)5(^-j3^), 



where A{x) and S{x) are the lump functions given by (2.4), and j denotes the 
j*^ lattice site located at fj = (ji, j2, is) £ ■ On account of (3.41), the V-A^ = 
condition gives 

a^^ - a], + h^^ - h]. + - c\,„ = (4.15) 
where the subscripts refer to the sites 

j = (ii,i2,i3), j- = (ii - l,i2, js) , 

j- = (31,32-1,33) and j'l' = (ji, ja, - 1) • (4.16) 

Equation (4.15) is the familiar field- flux conservation relation on a lattice (anal- 
ogous to the Kirchhoff law). In our case, this discrete form co-exists with its 
continuum realization V • A = at all f. 

(A comment on notation: Whenever the letter £ appears as a superscript, 
it denotes the color index. Otherwise, it may designate the lattice spacing, as in 
X — ji£ above, or it may indicate quantities associated with a lattice; e.g., the 
lattice coupling constant g£ that will be introduced later.) 

Figure 2 gives a simple example of a localized configuration with 
Ai = S{x)[A{y)-A{y-£)]A{z)S'', 



y 

and 



M = -[A{x)-A{x-£)]S{y)A{z)S'\ (4.17) 



Ai = 
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(the factor 5^^ can be replaced by any other distribution in the SU (n) index 
£). As explained in the caption, by examining the circulation of link fluxes 
associated with the four links of the central plaquette in Figure 2, one sees that 
(4.15) is fulfilled. Through direct differentiation, this configuration also satisfies 
V ■ = everywhere; it can be readily verifed that there are altogether 2x9 
cubic cells with nonzero . The number inside each of the nine plaquettes (all 
located at z = 0) in Figure 2 denotes the corresponding integral J B]. dxdy 
over the plaquette. 

4.-2. Lattice Variables in the Coordinate Space 

In the Bloch wave-number space, the generalized coordinates qi{K) in the 
expansion (4.10) satisfy 

<liiK) = qii-K)^ (4.18) 

and the transversality condition (4.12). At a given superscript £ and for each 
pair ±K 7^ , we can construct four independent Hermitian variables from the 
two combinations 

Ql{K) ± QliKy (4.19) 
each with the index t = 1 or 2, where 

QliK) ^ Y.e,{K),qi{K) (4.20) 

and et{K) given by (3.47). Thus, for a finite lattice of A/" = unit cubic cells 
with periodic boundary conditions, there are 2{M — 1) independent Hermitian 
variables for X ^ , and for K — Q 

3 constant solutions, (4-21) 

making a total of 

27^+1 (4.22) 

independent Hermitian variables in the Coulomb gauge, all explicitly exhibited. 
Here, as well as throughout this and the next subsections, we keep the group 
superscript £ fixed; therefore in the counting of degrees of freedom we do not 
include the trivial factor — 1 due to the SU{n) generators. 
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However, the same problem appears in a different form, if one wishes to 
use the coordinate space lattice variables. We note that the expansion (4.14) 
employs 3 jV Hermitian variables 



{a^,6j,c^.}. (4.23) 

Because of (4.15), there are M —1 constraints (in which the subtraction of 1 is 
on account of the sum of (4.15) over all j being identically); this leads to the 
same total number 3 jV — (jV — 1) = 2 jV + 1 independent variables, given by 
(4.22). 

Excluding the three constant K = Q solutions (mentioned in (4.21)), we 
may write 

= V X (4.24) 
where the spatial components of are given by 

ii = ^3 ^(^ - ^'i^) ^(y - ^2^) - ^3^) ' 
j 

4 = lYl S{x - h£) S{y - h£) A{z - h£) , 

j 

with the factor £ denoting the lattice spacing, as before. Clearly, and, 
therefore, also (4.23) are invariant under 

ji _^ /^ + y^^_ (426) 
For a lattice interpretation of these continuum relations, we specify 

= E x' -^(^ - -^'i^) ^(y - 3^^) - -^3^) • (4.27) 

Because of (3.41), (4.24) can be expressed in terms of the discrete variables: 



(4.28) 





= d 


- 


- rfj + ^/j. 










4 


= ^1 


-4. 
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similarly, the transformation (4.26) takes the discrete form 

~^ + Xj ~ Xj'_ ; 

V'j ^ V'j + Xj - Xj'^ , (4.29) 

Cj - Cj+Xj-Xj'^, 

where the lattice sites j, j'l and j'H are given by (4.16). 

As remarked before, we may identify aj, bj and as "link-fluxes" which 
flow along the three links connecting the site j to sites j', j" and j'" , given 
by (3.18). Consider the cubic cell in which x — ji£, y — and z — j^i are 
all between and £ . Label its plaquette shared by links jj" and jj" as Xj ; 
likewise the plaquette Yj is shared by links jj'" and jj' , and the plaquette 
Zj by jj' and jj" , as shown in Figure 3(a). Identify , ryj and (j as 
the circulating "plaquette-fluxes" which flow counterclockwise along the edges 
of plaquettes Xj ,Yj and Zj , as illustrated by Figure 3(b). (For example, 
resembles the current loop represented by a magnetic moment pointing in the x- 
direction, that is normal to Xj . The equation A — V xl resembles J — V x M 
from electromagnetism.) 

The continuum relation = V x implies V ■ A^ = . The discrete 
expression (4.28) carries the same implication: Each link is shared by four pla- 
quettes; through superposition of these plaquette-fluxes, the net flux along each 
link is given by (4.28), which automatically satisfles Kirchhoff's law (4.15). To 
summarize: in accordance with (4.21)-(4.22), there are 2J\f + 1 independent 
link- flux variables {a^bpCj}; excluding the three constant solutions, the 
remaining 2Af -|- 1 — 3 = 2J\f — 2 independent link- flux variables can all be gen- 
erated by the plaquette-flirxes through (4.28). In other words, among the SN" 
plaquette-flux variables { ryj, Cj} ■> there are 

3J\f-(2X-2) = U+2 (4.30) 

redundant ones that imply no net link-flux. (See Section 4.3 below.) The prob- 
lem of eliminating the redundant variables can then be shifted from link-fliixes 
to plaquette-fliixes. 

In the lattice formulation of QCD presented in this section, there is a precise 
connection relating the continuum equations (4.1), (4.24) and (4.26) with their 
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discrete counterparts given by (4.15), (4.28) and (4.29). The discrete reahzations 
can be further developed through the interplay between the lattice and its dual 
structure as follows: 

In the dual lattice, let links Xj, Pj, zj be duals to plaquettes Xj, Yj, Zj 
in the original lattice. The plaquette- fluxes ^j, ijp and C| can also be regarded 
as "currents" in the dual lattice, all flowing along the respective links. In the 
same dual lattice, assign to each lattice site j a "potential" x^j ; transformation 
(4.29) generates a corresponding current change along each link by an amount 
equal to the potential difference between its two end points. Such a change does 
not alter the link-fluxes (a^, Cj) in the original lattice. These extra degrees 
of freedom enable us to eliminate A/" — 1 of the N' + 2 redundant variables 
referred to in (4.30); the three remaining redundant degrees of freedom may be 
represented by adding arbitrary constants to ^j, ryj and i as we shall see. 

4-3. Elimination of Redundant Plaquette- Fltix Variables 

The lattice variables in the Bloch wave number space can be readily ex- 
pressed in terms of the 2J\f + 1 physical momentum-related ones, given by 
(4.20)- (4.22). However, for applications to confined configurations, such as glue- 
balls, etc., the use of lattice coordinate space variables gives a more direct physi- 
cal description. In this section, we give an explicit procedure for the elimination 
of j\/'-|-2 redundant variables among the SjV plaquette- fluxes r]j and , in 
accordance with (4.30). Again, the group superscript £ is kept fixed throughout 
this section. 

In the dual lattice construct a "tree" called T, by connecting arbitrarily 
J\f —1 links (without any closed cycles). Start from one of the end lattice sites 
of T and label it 0. The tree connects to another lattice site, say 1, through 
one of its links. Continue this way to move along the tree through all its links. 
Since each new link leads to a new lattice site, the J\f — 1 links of T must lead 
from the site to M —1 other lattice sites; i.e., the tree T connects all M sites 
in the dual lattice. 

Embed the dual lattice L d within the original cubic lattice L ; each of the 
links in Lo penetrates through a plaquette in L . Assign arbitrarily a distri- 
bution of 3J\f plaquette- fluxes (^|, ryj, and ) in each corresponding to a 
link-current in Lo . Again, start from the end site of T in ; by assigning 
the potential difference Xo ~ Xi between the sites and 1 to be the negative of 
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the link current along 1 , we can change that link-current (and therefore also 
the corresponding plaquette-flux in L) to zero through transformation (4.29). 
Proceed the same way along the tree T ; since there are M —1 arbitrary poten- 
tial differences between its M lattice sites, we can transform all the link-currents 
along T to become zero m. Ld', i.e., all the corresponding 

A/"-! (4.31) 

plaquette-fluxes in L can be set to zero through (4.29), without affecting any of 
the link- fluxes (aj, 6^, c^) . Since the lattice sites vn Ld are all located at the 
centers of the cubic cells in L , the entire tree T does not touch the surfaces of 
the original lattice L . In other words, none of the jV — 1 plaquette-fluxes that 
have been set to zero lie on the surfaces of L . 

Without the boundary conditions, the entire cubic lattice L would have 
six square surfaces. Take any one of these surfaces, say S ; it would have four 
edges. Assign an overall plaquette-flux flowing along the edges of S ; be- 
cause of periodic boundary conditions, we see that the edges facing each other 
become one, but with flowing along opposite directions, therefore cancelling 
each other. Furthermore, the six surfaces of L reduce to three different ones. 
An overall plaquette-flux on any of these three surfaces produces no link-fluxes. 
Thus, from the tree T, we can add three links in the dual lattice Ld , with 
each of the three different surfaces of L penetrated once: Assign zero also to 
the link-currents along these three additional links in Ld , therefore making the 
three corresponding plaquette-fluxes in L also zero. (This is a specific way of 
eliminating the three remaining redundant degrees of freedom by adding con- 
stants to ^j, Tjj and , mentioned at the end of Section 4.2.) Together with 
(4.31), we have succeeded in setting A/"— 1 + 3 = A/'+2 plaquette-fluxes to zero, 
without affecting the link-flux distributions in L . This then removes the re- 
dundancy mentioned in (4.30). The remaining 2M —2 nonzero plaquette-fiuxes 
(Cj; ^j; Cj) and the three constant of (4.21) give the 2J\f + 1 independent 
lattice coordinate space variables. 

An alternative procedure to eliminate the redundant lattice variables in the 
coordinate space is given in Appendix B. 
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5. LATTICE COUPLING CONSTANT 



5.1. Lattice QCD Hamiltonian 

For the continuum QCD, the Hamiltonian Ti , (4.3)-(4.4) in the Coulomb 
gauge, is written in terms of A^{f,t) and its conjugate momentum V^{f,t) , 
which satisfy (4.1) and (4.2). In our lattice formulation, A^{f,t) and V^{f^t) 
are replaced by A^(r, t) and n^(f, t) , given by (4.10)-(4.11), which denote the 
same field operators, but restricted to the O^'^-band Bloch functions. At a given 
lattice size I and for a finite total volume O = {N£)^ , the system has only a 
finite number of variables, as discussed in the previous section. Correspondingly, 
the continuum Hamiltonian 7i reduces to the lattice Hamiltonian 



^ J [J-'uiJUi + BiBi] d\ + Hconi (5.1) 
where, in place of (4.4)- (4.8), the lattice Coulomb interaction is 



^ 2 



^^coui = ^^|J-^y"cTi(r)(£,r|(V„D„)-^(-V2)(VbD6)-iK',f') 



Jaiif')d\d^r', 



(5.2) 



with the charge density now being 

4 = f'^^'A^^U:, (5.3) 
the color magnetic field given by 

eabc Bi = Va Ai - Vb Ai + g, f^^ A^ Al , (5.4) 
Da denoting the covariant derivative matrix whose elements are 

and J the O^'^-band Faddeev-Popov determinant (Jacobian) 

J = det I Va £>a I • (5.6) 

As in (4.1)-(4.8), all superscripts £,m,n refer to the group indices, from 1 to 
— 1 for SU (n) and all repeated indices are summed over. (Notice that the 
subscript £ in g£ denotes "lattice".) 
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The relation between the continuum and the lattice Hamiltonians Ti. and 
H may also be summarized as follows: Regard the continuum QCD Hamiltonian 
as a functional ^ of , and the bare coupling • We write 

n = J^{A',V',go) (5.7) 

in accordance with (4.3)-(4.8). The lattice QCD Hamiltonian H denotes the 
same functional J-' , but with and go replaced by A", 11" and gi; i.e., 

H = .F(l«,n«,^,). (5.8) 

Since the continuum QCD action is noncompact, so is the lattice QCD action. 
Because H consists of only a finite degree of freedom, the lattice coupling 
constant g£ does not have to be renormalized, unlike (4.9) for the continuum 
case. In the following, we shall examine the relation between the lattice coupling 
g£ , the continuum bare coupling cind its renormalization. 



5.2. Bare and Renormalized Continuum Coupling Constants 

A convenient and often used definition of the renormalized coupling in the 
continuum case is to express it in terms of the interaction energy associated with 
two external static color charges e{ and 63 , positioned at a distance R apart: 
Replace the continuum charge density (4.5) by 

^'(r) = f'^''A^{f)V:{r) + aUn (5.9) 

where 

4Ar) = e[ 5%T- R,) + ei 6'{f- 4) , (5.10) 

with e{ and 62 located at ^1 and R2 , and both being (n^ — l)-dimensional 
vectors in the group space. Let £^12 denote the part of the energy that is 
proportional to the product times 62 . Write 

— * — * „ 

where R= | i?i — i?2 | . At a given R , the coefficient g^^ can serve as the square 
of the renormalized coupling. In perturbation series® (assuming << 1 ) 
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2 2 

9r = % 



+ 0{gl) 



(5.12) 



where A is the ultraviolet momentum cutoff and 7 = 0.5772157 is the Eu- 
ler constant. (The reason for separating out the constant 7 from InAR is 
connected with (5.64) below.) The inverse relation of (5.12) may be written as 



2 2 
9o = 9r 



1 + 



247r2 



(ZnAi? + 7) 



+ 0{gl) 



(5.13) 



Thus, at two distances R and R' , to the same perturbative order the cor- 
responding renormalized couplings gn and gui satisfy the familiar asymptotic 
freedom relation^ 



1 



1 



9r' 9r 



lln 
24^ 



InR/R'. 



(5.14) 



It follows then, for R' < R, g\, < ; in particular, when R' /R — > , g\, — > . 



5.3. Lattice Coupling gi As the Renormalized Continuum Coupling Constant 

Next, we consider the lattice Hamiltonian (5.1)-(5.2). Replace the lattice 
charge density (5.3) by 

criif) = f'^'^A^^inKin+aUn (5.15) 

with cTf^^(f ) given by the external charge density identical to (5.10), and then 
evaluate the same interaction energy £'12 that is proportional to e{ times 62 . 
Write the lattice result as 

where, unlike (5.11)-(5.12), {g]i)L is a finite function of the lattice coupling gj . 
Because H does not have the symmetry of the continuum, but only that of a 
cubic lattice, {gji)L depends also on Ri,R2 and the lattice size £. We may 

— * — * 

select any i?i and R2 , and set 

{9l)L = gl; (5.17) 
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i.e., we choose the lattice couphng g^, so that the interaction energy £^12 (-Ri, R2) 
calculated from the continuum Hamiltonian Ti, and that from its 0*'*-band ap- 

— * 

proximation H are exactly equal for a particular pair of position vectors Ri 

— * 

and i?2 ■ Because of the cubic lattice symmetry, this equality extends to all 
pairs of Ri , R2 which can be reached from the selected ones through a lattice 
translation and a cubic rotation. 

The precise meaning of this equality can be stated in an alternative way: 
Express the continuum coupling constant renormalization formula (4.9) in the 
form 

9q = ge + Sg; (5.18) 

i.e., instead of using (5.13), regard the lattice coupling gi as the renormalized 
coupling for the continuum case. Expand and in terms of the complete 
set of Bloch functions { F^[K | r ) } , given by (3.72)-(3.73). Decompose 

A!^ = A^ + 5A^ and = ii^ + Sp^ (5.19) 

where 5A^ and 5V^ consist of all terms that depend on n 7^ 0-band Bloch 
functions. Substituting (5.18)-(5.19) into the functional (5.7), we can express 
the continuum 7i as a sum 

n = H + Sn, (5.20) 

where 6Ti. contains all terms that depend on 6A^ , SV^ and 6g . The condition 
for determining 5g is: For e{ and 62 located at the selected pair of position 

— * — * 

vectors Ri and R2 , the interaction energy £^12 calculated from H equals that 
using H ; i.e., 5H gives no correction. At different ^1 and R^ there will be, 
in general, a correction due to SH . 

Because of asymptotic freedom, there is a wide range of R and lattice size 
£, for which the equality (5.17) can lead to ^| << 1 , in which case corrections 
due to STi can be computed perturbatively in powers of gf . In the present 
noncompact lattice formulation, there is no need to take the limit £ — > , since 
the continuum solution is independent of the lattice size £ . 

In the following, we assume that g^ is indeed small; this enables us to relate 
perturbatively gf with , defined by (5.11). It is convenient to write 
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2 2 

9r = 9t 



1 + 



lln 

247r2 



( T + 7 + A + 5 



+ 0(^1) 



(5.21) 



where, besides the Euler constant 7 , there is another numerical constant A and 
a function 

5 = 6{e,Ri,R2) 



which satisfies 



hm 5(e,Ri,R2) = 0. 



Accordingly, (5.18) becomes 



9o = 9£ 



lln 



-^gt{lnM-X-5)+0{gt) 



247r^ 



(5.22) 
(5.23) 

(5.24) 



In the rest of the paper, we shall evaluate the constant A and derive the formula 
for 5 . As we shall see. 



A = 3.3559, 



The expression for 5 is given in Appendix C (see (C.33)). 



(5.25) 



5.4.. Power Series Expansion of the Lattice Hamiltonian 

In the 0*'*-band approximation, the Hamiltonian (5.1)-(5.2) can be ex- 
panded into a power series of the lattice coupling constant : 



H = Ho + Hint 

where Hq is ge independent and 

i^int = 9eHi+ g1 H2 + g£H3 + gj H4 + 



(5.26) 



(5.27) 



The field A^(f ) and its conjugate momentum n^(r ) are given by (4.10) and 
(4.11). At any given time t, they can also be written in terms of the creation 
and annihilation operators: 



K,t \^2uJt{K) 



a[{K) et{K)a Fa{K \ f) + a^(K)t et(K)« F„(K | fy 



(5.28) 
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and 



a\{K) tt{K)a Fa{K I f ) - a[{K)^ et{K)a Fa{K \ f) 



K,t 



(5.29) 

where et{K) and u!t{K) , with t = 1,2 (or + , — ) are given by (3.47) and 
(3.56), al{K) and al{K)^ satisfy the usual commutation relation 



(5.30) 



As before, K and K' are Bloch wave number vectors within the Brillouin zone. 
The operators a|(i?) and a|(i?)1^ are related to qi{K) and pi{K) of (4.10) 
and (4.11) by 



1 



4{K) tt{K)a + 4{-K)^ tt{-K)a 



and 



(5.31) 



4{K) tt{K)^-a[{-K)Ut{-K)a 



In terms of these annihilation and creation operators, the zeroth-order Hamil- 
tonian in (5.27) becomes, apart from an additive constant which may be set to 
zero. 



Ho ^ 0Jt{K) 4{K)^ 4{K) . (5.32) 

K,t,l 

The eigenvalues of Hq are 

K,t,l 

where A^l ^ = 0,1,2,-- - is the eigenvalue of the occupation- number operator 

al{K)'^ al{K) . Let the corresponding eigenstates be | AT > , with | > denoting 
the ground state of Hq . Thus, we have 
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4{K)\0>= (5.34) 

and Eo{0) = . To 0{gj) , the shift in the ground state energy Eo{Ri, R2) due 
to i/int of (5.26)-(5.27) is given by the famihar perturbation formula 



< TV I ffint I > 



|2 



Eo{R,. ^2) = < I i^int I > - ' ' ' + 0{Hl,) ; (5.35) 



the interaction energy £'12 , defined by (5.14), refers to the part of Eo{Ri, R2) 
which is proportional to the product ef 62 and corresponds to the following 
difference: 



E12 = Eo{RuR2)- lim Eo{RuR2) (5.36) 

where = | i?2 — -Ri | , as before. 

In the following, E12 will be expanded in powers of i^int : 

E12 = E^^ + E^;^ + • • • (5.37) 

where e'^2 is first-order in H\-at and E^^" second order in i^int . As we shall 
see, £^^2 contains 0{gf) ^ 0{gj) and higher order terms; likewise e[2^ contains 
0{gj) and higher-order terms. In this paper, we shall calculate £^12 only up to 
0{9t) . 



5.5. First Order in H\^t 

To first order in i?int , (5.35) reduces to 



Eo{Ri,R2) =<0|i/int|0>, (5.38) 

in which only the part proportional to ef 62 gives e[2 . Substituting the power 
series expansion (5.27) into (5.38), we see that since all odd gr^-power terms 
carry odd numbers of gauge field operators, whose vacuum expectation values 
are zero, the result is an even function of g(, . Because E12 is proportional to 
ef 62 , for our purpose we can equate i^int = -f^Coui ) which is given by (5.2). Up 
to 0{gf) , we can also set the Jacobian J = 1 . In power series of gi , we expand 
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(5.39) 



where Aa denotes the (n^ — 1) x (n^ — 1) matrix whose element is related 

to the 0*'*-band gauge field A^{f) by 

{AaY^' = f"^^' A"^. (5.40) 

Hence, we find 

= 9'i"i^+9tS^^^ + 0{gl) (5.41) 

where 

= -3efeJ <0\{e,Ri\V-^VaAaV-^VbAbV-^\e',R2)\0> . (5.42) 
From (5.28) and (5.34), it follows that 

< \Ai{r) A^{f') I >= J] [2utiK)] et{K)a et{K)b Fa{K \ r)Fb{K \ r')* 

K,t 

(5.43) 

where, as before, Fa{K \ r) is the 0*'*-band Bloch function, given by (3.38). In 
(5.43), the repeated subscripts are not to be summed over. Using the familiar 
expressions 



(i?l|V-2Va|f) 



(27r)3 p2 ^ 



(f I \R2) = / ^ e^P'-(^-'-^^) (5.44) 



and 



(27r)3 p' 

fflV-^lfM- [ }_piq-{r-r') 

^"^^^ '"^^"7 (27r)3g2e 



with p, and q denoting the magnitudes of p, p' and q, and combining these 
with (5.42)-(5.43), we derive^ 

= 3nef 4 E / [2-.™-^ 6,(K)„e,(K), ^^^^^^ 

X Fa(# I f ) F6(K I r'Y d^r d^r' 
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where 



(t> = p- {r- Ri) + q- (f ' - f ) + p'{R2 - 



(5.46) 



and the summation extends over the polarization index t = 1, 2 , K within the 
BriUouin zone, p,p' and q over a// wave numbers, and the subscripts a,b = x,y 
and z . 

Because e~^^"^ times the Bloch wave function Fa{K\r) has the periodic- 
ity of the lattice, the integrand in (5.45) consists of a periodic function in r and 
f' multiplied by a phase factor e^^P+^^-ir-r'} _ j^ggQiye the phase into a sum of 
terms 

(p + K -{f-f') = {p-q + K)-f-{p' -q + K)-f" + {p' ■R2-p-Ri) . (5.47) 
The integration over r and r ' gives 

p — q + K = 27rmi/£ 
and (5.48) 

p' -q + K = 27rm2/^ 

where rhi and m2 are three-dimensional vectors whose components are all 
integers. Introduce 

k = ^if +p) , 

^cm = i(i?2 + ^l) (5.49) 

and 

R = i?2 — Ri • 
The last term in (5.47) can be written as 

. -* -* 27r , ,-»-»-» , , 

p' ■R2-p-Ri = — (m2 -rhi) ■ Rcm + k- R. (5.50) 

Consequently, £^'^^ is a periodic function in with the periodicity of the 

lattice cell. 



48 



By using (5.48)-(5.49), we can convert the sum over p,p' and q into an 

— * 

integration over k times a sum over rhi and m2 ; this leads to 

where I{k) is a dimensionless function, given by 

^(^)= E j ^^.f^.^AKAe-'^^^^-^^^^- (5.52) 

mi,m2 

with the /c-integration over the entire continuum momentum space, but the 
i?-integration only within the Brillouin zone as indicated by the subscript B , 



X 



— 7 



(A;+I2l^7r)„(fc-2i™7rj;, 



(5.53) 



(A; + ^^i^ 7r)2 (A; - 7^)2 



and 02C„(^| m) given by (3.67) and (3.40), which is both dimensionless and 
independent of O. From (5.52), one sees directly that I{k) and therefore S^'^^ 
are both periodic in i?cm • 

To evaluate the contribution of S^'^^ to the constant A in (5.21), we examine 
the limiting case £/R ^ and call the result A*^*^ . The difference between the 
limiting result and that for an arbitrary nonzero £/R gives the corresponding 
contribution to the remainder S{i, Ri, R2) . When ^ , the integral (5.51) 
is dominated by the region kl = 0{R~^) . Furthermore, because of the rapid 



oscillatory phase factor exp i ^ (mi 
consider the special case 



TO2) • Rem. in (5.52), we need only 



mi 



m2 



m ; 



(5.54) 



— * — * 

as a result, I{k) becomes independent of -Rem. Prom (5.53), one can readily 
verify that as k and K both — > , the sum (5.52) over m is controlled by the 
single term mi = m2 = , for which the limiting behavior of the corresponding 
/mim2 {^1 ^) is given by the continuum function fc : 
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MK,k) ^ MK,k) 



E 



2K{k + K)^ ^ 



A;2 



(5.55) 



where K and k refer to the magnitudes of K and k , as before, and the sum 
over the polarization index t yields 



ab 



(5.56) 



To facilitate the integration (5.52), we decompose I{k) into a sum of two 
terms: 



where 



I{k) = Ic{k)+I'{k) 



(5.57) 
(5.58) 



and, on account of (5.54), 



J2frn,rn{K,k)-fa{K,k) 



(5.59) 



When £/R — > , the integration Ic{k) can be evaluated analytically. As shown 
in Appendix C, the result is 



^cik) = J^i-lnke + Xc) 



(5.60) 



in which the constant Xr is 



Xc = - + Inn / dx 



ln{x^ + 2) 



71" Jo (x^ + l)Vx^ + 2 



1.67040 



(5.61) 



The second term I'{k) in (5.57) is well-behaved at small k; its value at ^ = 
can be written as 



/'(O) 



67r2 



A' 



(5.62) 
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with the constant A' given by 

X' = -0.1496. (5.63) 

By using 

/fl^k ^ — 1 

^ ^ ilnke)e^'^-^ = [(/ni?/£) + 7], (5.64) 

we derive, as — > , 

^^^^ ^ ( 2^) ^ ^/^^ + ^ + ^^^^ ^ ^^-^^^ 

where 7 is the Euler constant and 

A« = Ac + A' = 1.5208. (5.66) 



5.6. Second Order in Hi^t 

In accordance with (5.35)- (5.37), 



I <0|iJint|iV> P , ^ I <0|gint|iV> P 

=-2l E^) ^tL2Z ' (5-67) 



in which the relevant states | N > are those of two gauge quanta, i.e. | N >= 
I 2 > with 

|2>= aJ^n^i)^«S(^2)+|0> . (5.68) 
To 0{gj) , we can write, similar to (5.41), 

E^^ = gt£^''^+0{gl) (5.69) 

and set, in (5.67), i^int = -f^Coui with J = 1 and = Va , so that the relevant 
part of HconX is 

i^Coui = g^i J d'r J d'r ' aiir) aUf') (f | - | r') (5.70) 
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where cr\{r) = /^^" ^^(f ) n^(f ) is the first term in (5.15). By using (5.28)- 
(5.29), (5.68)-(5.69) and the last equation in (5.44), we obtain 

< 2 I ifcoui \0>^ f'^'^ Yl ^^1*^ (^'^ ^2) {e{ e'^^'^^ + e"^^"^^) 

(5.71) 

where and £2,^2 denote the group and polarization indices of the two- 

gauge quanta in | 2 > , q \s the momentum variable in (5.44) for , which 
is related to the Bloch wave numbers Ki and K2 of the two quanta by 

q = #i+K2 + ^m, (5.72) 
with the components of rh all integers. The factor Mt-^t^ is given by 



2^ uJt,{Ki) LJt^{K2) ffl (5.73) 
X J d^rFa{Ki\ryFa{K2\n*e"^". 

We observe that when ti = ^2 = ^ , as Ki + K2 ^ 

Mt^tMKuK2) ^ (5.74) 
for all m including 0, since in this case 0012(^2)— LUtiiKi) = ijOt{K2)—ijOt{Ki) — > 

— * — * 

. Also, when ti ^ t2 but m = , as + K2 — > we have again 

Mt^t2{q;KuK2) ^ 0, (5.75) 

because 

J d^rFa{K^\rYFa{K2\rY ^ j d\ Fa{K^\rr Fa{K^\r) = 1 
which makes the subsequent sum over the polarization vectors in (5.73) 

o 

Both properties (5.74) and (5.75) will be useful for calculating £^(") in the limit 

52 



Substituting (5.71)-(5.73) into (5.67) and (5.69), we derive 



where 

(t>' = q-Ri-q'-R2, 

,-=i?. + i?. + H!L!5, (5.77) 

and the sum is over the polarization indices ti and t2 , the two Bloch wave 

number vectors Ki and K2 , as well as the vectors m and m ' , both having 
integer components. From (5.49) and (5.77), we have 

0' = {fh-m')-R^^ + k-R (5.78) 



where 



k = Ki + K2 + j{m + m'). (5.79) 



Hence, £^(") is also a periodic function in , the same as E^'^^ . Likewise, as 
£ — > , we can set 

rh = m' (5.80) 

and £^("^ becomes independent of • 

To derive the asymptotic behavior of in the limit £/R ^ , we follow 
steps parallel to those for . Write 



/jSj. „ik-R 
V ''^^^ ^'-''^ 



J I 

where 



II{k) = / J^HK^k) (5.82) 



B 



with the k integration extending over the entire momentum space but the K 
integration only within the Brillouin zone. The function h is given by 

Mi?, k)^'y: M..(g;K,KrM., (.iK,KO 3 
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where K' is a dependent variable, given by K' = {k — -K^)mod27r/£ ! i-^-, like 
i?, K' also lies within the Brillouin zone, but with the difference [K' — (k — 
K)]a — {2n/£)x integer. Here, as before, the subscript a = x,y and z. We 
note that as k ^ , the limit h{K, k) is well behaved, on account of (5.74)- 
(5.75). However, after the K integration, the function II{k) diverges as ^ — > . 

As in (5.57), write 

II{k) = IIc{k)+n'{k) (5.84) 

with 

im = ^ 1^ he{K, k) (5.85) 

and 

11' {k) = / [h{K,k)-K{K,k)]. (5.86) 

Jb {^'^r 

Choose hc{K, k) to be the integrand of II(k) without the O^'^-band approxi- 
mation, i.e. 



^k^ K\k-K\{K+\k-K\) 



(5.87) 



— * 

with the constraint that K lie within the Brillouin zone imposed as a boundary- 
condition on the integral (5.86). One can readily verify that II'{k) is free from 
divergence and obtain (see Appendix C for details) 

where the constant 

A(") = 11 in 2 + Ac + A" (5.89) 

with Ac given by (5.60) and 

X" = 247r^//'(0) = / 3 lim[/i(K,A;)-/ic(i^,A;)]. (5.90) 
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Note that in the hmit when the magnitude A; — > , the integrand remains 

A — * 

dependent on the direction k = k/k. (See (C.17), (C.25) and (C.28) below.) In 
Appendix C, we show that A" = 0.1238 and therefore 

A^^^) = 1.8362. (5.91) 
Substituting (5.88) into (5.81) and carrying out the integration over k , we find 

£:(^^) = (^/nf + 7 + A(-)) . (5.92) 

Combining (5.65)-(5.66) and (5.91)-(5.92), we estabhsh the constant A in (5.21) 
to be the one given by (5.25); i.e., 

A = A(*) + A^"^ = 3.3559 . (5.93) 
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APPENDIX A 



In the following we wish to compare the lattice formulation given in this 
paper to the conventional ones in which the lattice structure consists of only- 
discrete sites connected by links, but without being embedded within a contin- 
uum. An example of such a comparison for the case of a scalar field in one 
dimension was given at the end of Subsection 2.3. In this appendix we extend 
the comparison to the gauge theory. 

Let it be agreed that the aim is to discretize only the spatial coordinates 
while time remains continuous. This violates relativistic symmetries, but no 
fixed spatial lattice can be invariant under the full Lorentz group anyway. 

Take a simple cubic lattice of spacing I and period Nl in each orthogonal 
direction; the volume associated with the full period is O = {N£)^ . Conven- 
tionally, one labels the sites of the lattice as j = {ji, 32133) and attaches to 
each link a variable analogous to the vector potential, say a^, 6j, Cj on the links 
leaving site 3 in the a;— , y—, 2— directions (as in (3.12)). 

Consider first an Abelian gauge theory. The conventional lattice approach 
would be to write a lattice Lagrangian for a discrete system 



Cd = — {ctj — dji djY + ipj — djii + djY + (cj — dji'i + djY 



^ [i^j + bj' - - bjf + (bj + Cj" - hju, - Cjf + {cj + ay» - Cjf + 

(Al) 



a, 



where, as in (3.18), 



/ = (ji + I,j2,i3) , 

/' = (ii,i2 + 1,^3) , 
/" = (ii,i2,i3 + 1) ; 



(A2) 



the first term of (A.l) is "electric" and the second "magnetic". The new vari- 
able dj attached to site j is analogous to the "time" component of the vector 
potential. 

One notes that Cd is invariant under the lattice gauge transformation 
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CLj — ^ Clj -\- \ji — Aj , 

^ ^ ^ ^ (A.3) 
Cj — > Cj -\- Xj'" — Xj 5 

dj dj + Xj 

where Xj is an arbitrary quantity attached to site j . One may take advantage 
of this to pass either to the lattice time-axial gauge where 

dj = (aU j) (A.4) 

or to the lattice Coulomb gauge where 

aj - aj'_ + bj - bj'^ + Cj - Cj>>> = ; (A.5) 

here, as in (4.16), 

j- = Ui - 1,J2,J3) , 

j'l = (ii,j2-l,j3), (A.6) 
j- = (Ji + l,j2,j3-l)- 

In time-axial gauge one may write the Lagrangian (A.l) in "momentum 
space" : 

K 

where q = q{K) and u = u{K) are 3x1 column matrices. The components 
of q are the Fourier coefficients of a^, 6^, Cj . Those of u are 

MO) = ^ si^Y' (A-8) 

— * — * 

with 9 = K£ . Evidently there are redundant degrees of freedom, which one 
may eliminate by writing 

Qi{K) = u{K) ■ q{K) (A.9) 

where (u)a = Ua / "n/u^, and 

Qt{K) = et{K) ■ q{K) {t = 1,2) (A.IO) 
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with u, ii, 62 forming a right-handed orthonormal system. Then (A. 7) becomes 

= Y1 (qI Qe + Yl (Ql Qt - Ql Qt)] ■ (A.ii) 

In this gauge there is a supplementary condition on the state vector, 

^ < Ql > = (A.12) 



dQi{K) 



in the absence of sources, derived from the lattice analog of Gauss' law. Hence 
the variables can be ignored and the Hamiltonian written as 

2 

K *=1 

where Pt{K) is the conjugate momentum of Qt{K) . 

Alternatively, in lattice Coulomb gauge one has from the outset = , all 
K ^ Q ^ and the variables dj are decoupled from Qi . In the absence of sources 
the dj vanish and the Hamiltonian again takes the form (A. 13). 

All this proceeds smoothly, but one encounters difficulties in extending it 
to a non-Abelian theory. For example, how is one to represent the nonlinear 
term in the magnetic field? 

Two different approaches are possible. The ffist is to go over to a compact 
formulation^ based not on aj, 6j, Cj but on e*^"-' , e*^''^ , e^^^^ . The disadvantage 
is that this formulation, for nonzero £ , is no longer connected with the true 
continuum theory based on 

where the are continuous fields with gauge transformation 

A,^A, + ^. (A.15) 

The "compact" theory derived by nonlinearizing (A.l) and (A. 3) is an abstract 
model of the original theory, not a concrete approximation to it. 

The second approach, which we pursue in this paper, is to seek a noncom- 
pact lattice theory that can serve as a first approximation in a systematic way 
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which can approach the exact continuum theory, without requiring I to ap- 
proach zero. Thus (A. 14) and (A. 15) will be our starting point, not (A.l) and 
(A. 3); the lattice variables aj,hj,Cj will enter only as parameters controlling 
the few degrees of freedom allowed to A in the zeroth approximation. 

— * 

A natural way to implement this idea is to make each component of ^ a 
constant along any link in the corresponding direction, and interpolate linearly 
between links. Referring to the lump functions defined above, of which Li is a 
constant and L2 is linear, we see that this amounts to setting 

A{r) = ^[a,xLi(x-ji£)L2(y-j2^)L2(^-j3^) 



+ bj yL2{x- jii) Li{y- ^2^) 1^2 - jsi) 
+ Cj zL2{x- jit) L2{y - j2i) Li{z - jd) ] ■ 



(A.16) 



The idea is to substitute this expression into (A. 14) so as to obtain a Lagrangian 
depending on the "effective" lattice variables aj^bj^cj. (Writing A(f),Li{x) 
and L2{x) as V{f) , C{x) and A{x) , we see that (A.16) is the same as (3.7)). 

For an Abelian theory this idea works nicely in time-axial gauge because 
the time-independent gauge transformation (A. 3) actually generates a transfor- 
mation of the form (A. 15) on the continuous fields defined in (A.16). As shown 
in (3.22), this makes it possible to implement the subsidiary state condition 
derived from the continuum Gauss' law by eliminating some of the lattice vari- 
ables similar to (A. 13). But this no longer works in a non- Abelian theory: the 
non- Abelian version of (A.l) is not invariant under any transformation of the 
type (A.3). 

For the non- Abelian theory, we propose instead to work in Coulomb gauge 
(or Coulomb-like gauge (3.84)). We find that if we replace the function Li 
in (A.16) by L3 , then on account of (2.3) the condition (A. 5) on the lattice 
variables is sufficient to ensure the true Coulomb gauge condition 

V-A = (A.17) 

at every point in space. Moreover, if the higher band functions are built in the 
manner indicated by (3.73), then (A.17) can be satisfied independently within 
each band at every point in space. This remains true in non-Abelian theory 
because the Coulomb gauge condition is linear. The details are carried out in 



59 



Subsection 3.2 and Section 4. Identical considerations can be applied to any 
gauge in which the gauge- fixing condition has no nonlinear terms, as in (3.84). 
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APPENDIX B 



In this appendix, we give another method for the removal of redundant 
lattice coordinate space variables, alternative to the one given in Section 4.3. 

The objective is to eliminate the M -\-2 redundant variables among the set 



U'i.rilC^] (B.l) 

of 3A/^ plaquette flux variables. Again, the group superscript I will be kept 
fixed, and the factor of its ■n? — 1 degrees of freedom due to SU (n) will not be 
included. 



Theorem. For any plaquette flux distribution (B.l), through (4.29) in which 

{e.^vlCU - {ef,<,Cf} (B.2) 

the transformed plaquette fluxes can be made to satisfy the lattice divergence- 
free equation: 

C^-ef + ry^'-ryf + C^t-Cf = (B.3) 

at all j , with the sites j' , j" and j'" given by (3.18). 



Proof. Consider the cubic cell Tj that has Aj, Yj and Zj of Figure 3(a) 
as three of its six surfaces. The other three surfaces are (in the same notation) 
Aj/, Yjii and Zjm . Define the lattice differences 

dy (lyh ^ 4' - , (B.4) 

and 

to be the net flux flowing out from Tj along the x^y and z axes. The lattice 
divergence 

a-// ^ Y.da{li), (B.5) 
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represents the total flirx leaving tj . 

Next, we turn to the dual lattice Ld The center of Tj in the original 
lattice L becomes the lattice site j in Ld ■ Correspondingly, ij^ and 
denote link currents, and d ■ Ij^ is the net current leaving the site j (with 
the components of Tj^ now defined on the links j'j, j"j and in Ld). 
Likewise, introduce 

d^x- = Xj"-2x-+Xf:. (B.6) 
Oz Xj = Xj'" ^Xj + Xj'^' 

and 

d'^^Y^dlx] (B.7) 

to be the second-order lattice differences and the lattice Laplacian. The sites 
J, j', ■ ■ j'l! m Ld are dual to the cubic cells Tj, tj/, ■ ■ tj'" in L which, in turn, 
are related to the sites j, j', ■ ■ j'H in L , given by (3.18) and (4.16), in the way 
described above. 

Equation (B.3) states that through the transformation (4.29), Ij ^ becomes 
I'j- which satisfies d ■ I'j^ — ; therefore 

-d' ^ d-l/ ^ . (B.8) 
We can readily obtain the solution through the lattice Green's function 

where Qa = -f^a^ as before, and the sum extends over the Brillouin zone but 
with its origin excluded. The result is 

= E '5{?i-n)p\. (B.io) 

i 

The theorem is then established. 
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The function Q{f) has been well studied in the literature. For example, at 
r = and — > oo 

^(0) = 0.25273, (B.ll) 
which can also be expressed^'^ in terms of the elliptic integral 

^(0) = ^ (18 + 12V2-10V3-7v^) X K'^[{2 - Vs) (3 - V2)] (B.12) 
where 

K{k) = [\l-k^sm^(t))'^ dcf). (B.13) 
Jo 

Returning to the problem of eliminating redundant variables, we note that 
(B.3) imposes M — 1 conditions. Since any plaquette flux distribution can be 
transformed into one that satisfies (B.3), the resulting divergence-free distribu- 
tion consists only of 

3J\f-{H-l) = 2J^+l (B.14) 

independent plaquette-flux variables. As explained at the end of Section 4.3, 
we can arbitrarily set three surface plaquette fluxes of L to zero (one along 
each Cartesian direction) without affecting the link- flux distribution ( a^, bj, Cj). 
This and (B.14) complete the reduction process: the number of independent 
variables in (^j, ryj, (j) becomes 2jV -|- 1 — 3 = 2Af — 2 , in accordance with 
(4.30). 
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APPENDIX C 



C.l. The Integral Ic{k) 

To calculate the integral Ic{k) given by (5.58) for kl « 1 , we divide the 
Brillouin zone, — f- < Ka < j (where a = x,y,z) , into two regions: the interior 
of a small sphere s with its center at the origin K — Q and the volume outside 
the sphere. The radius of the sphere, Kq , is chosen such that k « Kq « j . 
Correspondingly, the integral (5.58) can be separated into two terms: 



Ic{k) = lin + l 



out 



where 



and 



[ ^ MK, k) 

J in 



4ut = / JTT-^ fc{K, k) 
out K^T^) 



(C.l) 
(C.2) 

(C.3) 



with fc{K,k) given by (5.55), the integral (C.2) referring to the integration 
inside the sphere s and (C.3) for the integration outside s but within the Bril- 
louin zone. Since /in contains Inki , k may not be set to zero in its integrand. 
A straightforward calculation in spherical coordinates yields: 



(C.4) 



where we have dropped the terms that vanish in the limit A; — > . To the same 
order, k may be set to zero in the integrand of /out , on account of the low 
momentum cutoff {K > Kq) . We have 



lout — 



'out (27r)3 fc-o 
It follows from (5.55) and (5.56) that 



lim^ fc{K,k). 



(C.5) 



lim/c(i^,^) = ^ 



1 - 



{k-Kf 



(C.6) 



Because of cubic symmetry, we may replace {k-K)"^ / K"^ by | in the integration; 
(C.5) becomes 
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1 f (f 



'out 



K 



3 7,,, (271)3 /.3 



(C.7) 





in which we introduce the dimensionless vector k = - K , whose Cartesian coor- 
dinates will be denoted by {x, y, z) , so that its magnitude is k = ^/ x'^ + y'^ + z'^ . 
The integration domain of (C.7) in terms of k is given by —l<x,y,z<l and 
n > Kq . In terms of spherical coordinates with a and b denoting the polar 
and azimuthal angles, we obtain 



lout 



247r3 



/•2TT PIT £ 

/ dh I da sin a In K,{a,b) ^ In (C.8) 

Jo Jo DTT^ TT 



where K{a, b) is the radial coordinate of a point on the boundary of the Brillouin 
zone in the direction specified by the angles a and b . The Brillouin zone is a 
cube bounded by six faces: x = ±1 , y = ±1 and z = ±1 . Because of cubic 
symmetry, we need only integrate over one of these faces and then multiply the 
result by six. Take the surface at z = 1: we have 



X = tan a cos b , 
y = tana sin 6, 



(C.9) 



and 



(C.8) reduces to 



K{a,b) = sec a; (C.IO) 



/out = ^ / dx dy -\ ^ + ln—- C.ll 

Carrying out the integration over y , we find 



Iout = -^[ln-— + l-- dx -—A-—^== . C.12 
Combining (C.4) and (C.12), we establish (5.60) and (5.61). 
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C.2. The Integral IIc{k) 

Following the same method as in the above section, we separate the integral 
(5.85) into two terms, 

Ilcik) = //in + //out (C.13) 

where, as in (C.1)-(C.3), 

Ilin = [ ^h,{K,k) (C.14) 



in 



is the integration inside the same small sphere s and 

//out = / -j:^hc{K,k) (C.15) 

^out {■^'^) 

is the integration outside s but inside the Brillouin zone. The integrand hc{K,k) 
is given by (5.87). In the limit k ^ we have 

//out = / ^ lim hc{K,k) = i /out (C.16) 

Jout l^'^j '^^O ^ 

since, on account of (5.87), 

lim MK,fc) = 1 (C.17) 
which can be replaced by {12K^)~^ because of cubic symmetry. Therefore 



//out = TTT^ In -— + 1 / dx ^ , ' (C.18) 



where, as before, /Cq is the radius of the small sphere s . 

To carry out the integration inside s for //in , we find it convenient to 
introduce the following prolate spheroidal coordinates for the Bloch wave number 

— * 

vector K : 

= ^ V(e'-l)(l-r/2) COS0, 
k 



Ky = - Vie-m-V^) sin0, (C.19) 
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with O<0<27r, — 1<?7<1 and ^ < 1 . To derive the In Ko/k term and the 
constant term of //in in the hmit /c , we may approximate the spherical 
domain of the integration by a prolate spheroid with foci (0,0,0) and ( 0, 0, /c ) 
and with the length of its semi-major axis equal to Kq , which corresponds to the 
upper limit of ^ < in the ^-integration. The error due to the approximation 
vanishes in the limit /c — > . We find 



IIw 



2Ko/k 



327r 
1 

24^ 



dr] — 
-1 ? 



k 



1 + 



-9^2- 



(C.20) 



Combining (C.18) and (C.20), we derive the expression of IIc{k) when k£ ^ 0; 
this leads to (5.88) and (5.89), except for the constant A" which will be evaluated 
below. 



C.3. The constants X' and X" 

The expression for the constant A' follows from (5.53), (5.54), (5.55), (5.59) 
and (5.62) and it reads: 



A' = 27r' 



d^K 

(27r)3 



E 

a.t 



nCa{9\m) 



K + 



27rm\ 

e ) 



1 



. (C.21) 



To derive the expression for A" , we separate the summation ^^^^^ 
h{K, k) of (5.83) into a sum over the same polarization, ti = t2 = t , and 
another sum over different polarizations ti ^ t2 ■, i.e. 



h{K,k) = hi{K,k) +h2{K,k) (C.22) 

where 

h.iK, k)='y: Muik;K,KrM,,ik;K,K') 
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and 



1 V- Mt,tAk;K,Kr Mt,t,ik;K,K') 



(C.24) 



with K' defined by the prescription following (5.83). The limit /c ^ of 
hi{K, k) can be readily obtained, and we have 



lim = iE'*'^""''' 



4 V ^KK) 

In the same limit, using the Fourier expansion (3.66) and (3.67), we find 



(C.25) 



j d^rFa{K\ryFa{K'\rye'^-^ = N^'^Ca{9\m)Ca{9-ke\m) = l + 0{k'^e'^) 

rh 

(C.26) 



and therefore 



^ ^ [u, (i?) - g;,, (K) ] ^ ^^^^ e., (-i?)„ + 0(fc^£^) 

(C.27) 

for ti^t2i in accordance with (5.75). Substituting (C.27) into (C.24) we obtain 



lim h2{K,k) = \ y 



[^,,(i?)-u;i,(i?) 



^tiW^t.(i^) [a;,,(iC)+a;,,(K)] 
X E [ eti (^)a ^ • et, (-K)„ ] [ et, (K)6 ^ et, {-K)b] . 

a,b 

(C.28) 



Following the definitions (5.90) and (5.86) we have 



A" = X'( + X'^ 



(C.29) 



where 



A'/ = 247r^ 



= 



(27r)3 fc^o 



lim [hi{K,k) - K{K,k) \ 



E 



^ROOt? 1 



2a;|(K) 



(C.30) 



68 



and 



a,b,c 



Both integrations (C.30) and (C.31) are well defined. The result is A" = 0.1238 , 



C.4- The Remainder 5 

The function S{£, Ri, R2) is defined by (5.21)-(5.23), which relates (7^ to 
g"^ . From (5.16)-(5.17), we see that appears as the coefficient of e\ 62 /47ri? 
in £^12 ; the latter can in turn be written as (neglecting 0{g^) ) 

£^12 = gf^+gtis^'^ + s^^'^). (C.32) 



It follows then 



S{£, R,, R2) = {S^^ + S^"^) -In^-^-X (C.33) 

1 i-Tl 6^ ^2 ^ 



where 7 is the Euler constant, A is given by (5.25), 

/j3i. pik-R _ 

with l(k) given by (5.52), and £^^") given by a similar expression (5.76). When 
^^0,5 satisfies 

lim (5(^,^1, ^2) = 0, (C.35) 

since that is the condition for the determination of A . In general, the dependence 
of 5 on i?i and R2 is not a simple one, though well defined through (C.33). 

Because of (5.50) and (5.78), both £^'^^ and are periodic in 

-Rem = 2 (-^1 -^2) 
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with the periodicity of the lattice. It is convenient to average over a unit 

lattice cell; this gives a < 5 >av which depends only on £ and R = R2 — Ri ■ 
For < 5 >av 5 we may set 

rhi = fhi and in = in' (C.36) 

in (5.52)-(5.53) and (5.76)-(5.77). Thus, the integrand in (C.34) can be written 
as 

^(^) = E / 7^fr^rn{K,k) (C.37) 

where 

frnrn{K,k) = /c'^ 17^ J] i ^. (C.38) 

Likewise, we have 

n(k) (C.39) 

where 

im = J^§^ M^, k) (C.40) 

and h{K, k) is given by (5.83). Substituting the above £^'^^ and f^") into 
(C.33), we derive the expression for < 5 >av ■ 
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FIGURE CAPTIONS 



Figure 1. Equations (3.3)-(3.5) give an electromagnetic field configuration in 
the time-axial gauge within the unit cell in which y, z all lie between and 

Figure 2. A picture of the localized QCD gauge field configuration (4.17) in 
the Coulomb gauge. All plaquettes are located on the {x, y) plane at 2; = . 
Move along the arrow direction of the central plaquette: Start from its lower 
left lattice site and label it 0; then label consecutively the other three lattice 
sites 1, 2 and 3. The gluon field described by (4.17) corresponds to a link- flux 
distribution (in the notation of (4.14)) aQ = h\ = —a\ = —60 = 1 and all other 
aj , 6^ and zero. Thus, Kirchhoff's law (4.15) is satisfied. The number inside 
each plaquette denotes the integral of f Bjdxdy over the plaquette. On 
account of the factor A{z) , the gauge field distribution (4.17) is nonzero over 
two layers of 9 cubic cells each (making a total of 18 cells), extending from the 
9 plaquettes at 2 = in the figure to z — ±£ . 

Figure 3. The lattice sites j, j', /' and j'" are defined by (3.18). (a) shows 
the plaquettes Xj ,Yj and Zj that are perpendicular to the x,y and z axes, 
(b) gives the plaquette-fiuxes , rjj and associated with these plaquettes. 
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